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Abstract 

Some Kharitonov-like robust Hurwitz stability criteria are established for a class 
of complex polynomial families with nonlinearly correlated perturbations. These 
results are extended to the polynomial matrix case and non-interval D-stability 
case. Applications of these results in testing of robust strict positive realness of real 
and complex interval transfer function families are also presented. 
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1 Introduction 

Motivated by the seminal theorem of Kharitonov on robust stability of interval 
polynomials [111 a number of papers on robustness analysis of uncertain systems have 
been published in the past few years|| f|, ||, || |7|, ||, |9|, [UJ. Kharitonov's theorem states 
that the Hurwitz stability of the real (or complex) interval polynomial family can be guar- 
anteed by the Hurwitz stability of four (or eight) prescribed critical vertex polynomials in 
this family. This result is significant since it reduces checking stability of infinitely many 
polynomials to checking stability of finitely many polynomials, and the number of critical 
vertex polynomials need to be checked is independent of the order of the polynomial fam- 
ily. An important extension of Kharitonov's theorem is the edge theorem discovered by 
Bartlett, Hollot and Huang |4j]. The edge theorem states that the stability of a polytope of 
polynomials can be guaranteed by the stability of its one-dimensional exposed edge poly- 
nomials. The significance of the edge theorem is that it allows some (affine) dependency 
among polynomial coefficients, and applies to more general stability regions, e.g., unit 
circle, left sector, shifted half plane, hyperbola region, etc. When the dependency among 
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polynomial coefficients is nonlinear, however, Ackermann shows that checking a subset of 
a polynomial family generally can not guarantee the stability of the entire family JIT, [12 



In this paper, we consider a class of complex polynomial families with nonlinear coef- 
ficient dependency. Based on our previous results, we will establish some Kharitonov-like 
robust stability criteria, i.e., the entire family is stable if and only if some critical vertices 
in this family are stable, and the number of critical vertices is independent of the order 
of the polynomial family. We will then extend our results to the polynomial matrix case 
and non-interval D-stability case. Applications of these results in testing strict positive 
realness of interval transfer function family are also presented. 



2 Main Results 

A polynomial p(s) is said to be Hurwitz stable, denoted by p(s) G H, if all its roots lie 
within the open left half of the complex plane C. A polynomial family P is said to be 
Hurwitz stable, denoted by P C H, if all polynomials in P are Hurwitz stable. 
Consider the n-th order real interval polynomial family 



r = jp(s) | p(s) = ^ > ft e h > ^ + ] > i = 0, 1, • • • , n j (1) 
and define the four Kharitonov polynomials of T as 

= % + + 4s 2 + qts 3 + qIs A + q^s 5 + ■■■ (2) 

K 2 {s) = q+ + qts + q^s 2 + q^s 3 + qjs 4 + g 5 + s 5 + • • • (3) 

K 3 (s) = qt + QiS + Q2S 2 + qts 3 + qts A + q^s 5 + ■■■ (4) 

K 4 (s) = go - + qfs + qts 2 + q^s 3 + q^s 4 + qts 5 + ■■■ (5) 

Lemma l(Kharitonov's Theorem for Real Polynomials) [|I| 

rcff^ K x {s) , K 2 (s) , K 3 (s) , K 4 (s) G H (6) 
Consider the n-th order complex interval polynomial family 



&=\8{s)\8{a) = 5^(a< + j/fys* , on G [a { , a+] , A 6 [ft , A + ] , % = 0, 1, • • • , n 
I i=o 



(7) 



and define the eight Kharitonov polynomials of A as 

K+(s) = {a Q +j^)+{a^+jPt)sHat+3Pt)s 2 Hat+3^)s 3 H^+3M^H^ 



(8) 
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K+(s) = (a +jp+)+(af+jpf)s+(a++jfe)s 2 +(az+jpz)s 3 ^a4+jpt)s 4 Hat+jPt)s 5 +- ' ' 

(9) 

K+(s) = (a++j(3 )+(a^+j^)s+(a 2 +j(3+)s 2 +(a++j(3+)s 3 +(a++j^)s 4 +(a^+jp^)s 5 +- ■ ■ 

(10) 

K+(s) = (a++jp+)+(af+j^)s+(a 2 +j^)s 2 +(a^+jf3+)s 3 +(a++jP+)s 4 +(at+jP5)s 5 +- ■ ■ 

(11) 

K{(s) = (ao+j(3 )+(at+j^)s+(ai+j(3+)s 2 +(a^+j(3+)s 3 +(a^+j^)s 4 +(ai+jp^)s 5 +- ■ ■ 

(12) 

K 2 {s) = (a +jp+)+(ai+j^)s+(a++j^)s 2 +(ai+jPt)^+(^+jPt)s i +(^+jp5)s 5 +- ■ ■ 

(13) 

K^(s) = (a++jP )+(a++jP+)s+(a 2 +j(3+)s 2 +(a^+j(3^)s 3 +(a++j^)s i +(at+jPt)^+- ' ' 

(14) 

K±(s) = (a++jP+)+(a^+jPf)s+(a 2 +j^)s 2 +(a++j^)s 3 +(ai+jp+)s i +(a^+jp+)s 5 +- ■ ■ 

(15) 

Lemma 2 (Kharitonov' s Theorem for Complex Polynomials) 

AcH^ K+(s) , K+(s) , K+(s) , Kt(s) , K{(s) , K 2 {s) , K^(s) , K^(s) e H 

(16) 

Now consider the n u -th , n„-th order real interval polynomial families T u and r t , . De- 
note their Kharitonov polynomials as K^(s) , i = 1,2,3,4 and Kj(s) , j = 1,2,3,4 
respectively. 

Similarly, consider the n u -th , n^-th order complex interval polynomial families A u 
and A„ . Denote their Kharitonov polynomials as Kf u (s) , K^ u (s) , i = 1,2,3,4 and 
Kf°(s) , Kj v {s) , j = 1,2,3,4 respectively. 

For any function f(x,y), define 

) = {/(Pu(s) , Pv(s)) | p u (s) e T u , p v (s) e T v } (17) 
f(A u , A v ) = {f(6 u (s) , 8 v (s)) | 8 u (s) e A u , 8 v (s) e A v } (18) 

Lemma 30 
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For any fixed complex number z G C, suppose the polynomial family T u — zT v has a 
fixed order. Then 

T u - zT v C H 

(19) 

Kf(s)-zK](s)eH, i,j = 1,2,3,4 

If the location of z is known, then the number of critical vertices need to be checked 
can further be reduced. For example, if z is on the negative real axis, then only 4 out of 
the 16 critical vertices need to be checked, namely 

T u - zT v C H <t=^ 

(20) 

K?(s) - zK?(s) G H , i = l,2,3,4; 

if z is on the imaginary axis, then only 8 critical vertices need to be checked; if z is in the 
left half of the complex plane, then only 12 critical vertices need to be checked^], |5], [7], |^]. 

For complex polynomials, we have the following similar result 
Lemma 4 

For any fixed complex number z G C, suppose the polynomial family A u — zA v has a 
fixed order. Then 

A u - zA v c H -<=>- 

(21) 

Kt u {s)-zKp{s), Kr(s)-zKr«( s )eH, ij = 1,2,3,4 

Theorem 1 

Consider the polynomial family 

a m r- + a m ^- l r v + a m _ 2 r™- 2 r 2 v + + a 2 r 2 u r™- 2 + oirj;- 1 + a r™ (22) 

where G R , k — 0, 1, • • • , m . Suppose it has a fixed order. Then 

- 1 + a r™ Cl^ 

a m [Kf{s)} m + a m . 1 [Kt{s)] m - 1 K]{s) + (23) 

+a 1 K?(s)[Kv(s)} m - 1 + a [Kv(s)} m G H , i,j = 1,2,3,4 

Proof: Consider the polynomial 

q(z) = a m z m + a m _iz m_1 + a m _ 2 z m ~ 2 + + a 2 z 2 + a x z + a (24) 

Let r = max{fc | ^ 0}. Then q(z) can be expressed as 

q(z) = a r (z - zi)(z - z 2 ) (z - z r -\){z - z r ) (25) 
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where z± , z 2 , , z r -\ , z T G C. Hence, we have 



+ + cuTj:™- 1 + a r™ C H 



v 



a r + a r _i (^) r 1 + + ai (^) + a 



C if 



a - (fc - *0 (fe - *») (fr " ^r-l) (ft " 



C H 



^r™ r (r„ — 2ir„)(r u — z 2 v v ) (r u — ^ r _ 1 r 1 ,)(r u — z r r„) c if 



T„ - C if , k = 1, 2, , r - 1, r 



r = m 



T u — z k T v C H , k = 1, 2, , r — 1, r and T v C H 



r < m 



Lemmas 1&3 



/q*(s) - z k KV(s) 6 if , i, j = 1, 2, 3, 4 , fc = 1, 2, ■ ■ • , r - 1, r 



K^{s)-z k KJ{s) e if, i, j = 1,2,3,4, fc = 1,2, -l,r and ifj(s) 6 if , j = 1, 2, 3, 4 



• • • - - e h , i, j = 1, 2, 3, 4 



[KV(s)]> 



K±[s) 
K v As) 



— Z r ^i 



K?(s) 



+ a r 



Kf{s) 



r-l I KV(s) 



Z 2 ■■■ 



EH, i,j = 1,2,3,4 



r-l 



+ 



+ ai 



+ ao 



Gif, i,j = 1,2,3,4 



a m [Kf{s)} m + a ro _i[^(s)] ro - 1 iq'(s) + 



• • • + ^Kfis) [K](s)] m ~ 1 + a [KJ(s)] m i,j = 1,2,3,4 

This completes the proof. 



From the proof of Theorem 1 and by Lemmas 2 and 4, we have 
Theorem 2 



Consider the polynomial family 



c m A™ + c m _! 1 A„ + c m ^- 2 A 2 v + + c 2 A 2 u A™~ 2 + Cl A U A™- 1 + c A™ (28) 

where CfcGC,/c = 0,l,---,m. Suppose it has a fixed order. Then 
c m A™ + c m _! A™- 1 A, + + dA^A™- 1 + c A- c H 

c m [Kt u (s)] m + c m ^[Kt u {s)r- l Kl v {s) + 

+c 1 Kt u (s)[K^(s)} m - 1 + c [K^(s)} m e H, (29) 

c m [K- u {s)} m + c^K-^^Kj^s) + 

+c 1 K- u (s)[K- v (s)r- 1 + c [K- v (s)r e H, i,j = 1,2,3,4 

Remark. We have established strong Kharitonov-like criteria for the stability of a class 
of polynomial families with nonlinearly correlated perturbations. The number of critical 
polynomials need to be checked is independent of the order of the polynomial family. 
Example 1 

Consider a negative unity feedback system with the forward path as three same blocks 
in tandem. Each block consists of an interval plant with negative unity feedback. 
Then, the characteristic polynomial of the closed-loop system is 

[N(s)f + [N(s) + D(s)} s (30) 

By Theorem 1, we only need to check 16 vertex systems for the stability of the entire 
uncertain system family. Furthermore, since all the roots of 

q(z) = 2z 3 + 3z 2 + 3z + 1 (31) 

lie within the left half of the complex plane, only 12 out of the 16 vertex systems need to 
be checked to verify robust stability of the entire system family. 
Example 2 

Consider a negative unity feedback system with the forward path as a controller and 
an interval plant in tandem. The controller is simply a gain k , but can be switched 

among {k\ , k 2 , , k m } under different working conditions. Thus, robust stability of 

the entire system family is tantamount to 

[hN(s) + D(s)][k 2 N(s) + D(s)} [k m N(s) + D(s)} C H (32) 

By Theorem 1, we only need to check 16 vertex systems for the stability of the entire 
uncertain system family. Furthermore, since all the roots of 

q(z) = (hz+ l)(k 2 z+ 1) (k m z + l) (33) 
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lie on the real axis, only 8 out of the 16 vertex systems need to be checked. Moreover, if 

k\,k 2 , ,k m have the same sign, then only 4 out of the 16 vertex systems need to be 

checked. 



3 Some Extensions 

3.1 Extension to Non-Interval D- Stability Case 

Given any stability region D in the complex plane C, a polynomial p(s) is said to be 
D-stable, denoted by p(s) G D, if all its roots lie within D. A polynomial family P is said 
to be D-stable, denoted by P C D, if all polynomials in P are D-stable. 
Let the uncertainty bounding set (hyperbox) be 



(34) 



(35) 



Q = {<? = (?i > 12 , , qi) T | ft e [q { , qf] , 

« = 1,2, ,0 

and define its one-dimensional edge set as 

Qe = {<? = (<2i , 92 , , qi) T I Qk e [q^ , qf] for some 

k e {1,2, , 1} and q t g {q~ , qf} for all t ^ fc} 

Consider the ni-th, n2-th order complex polynomials 

ni 

n( S ,g) = 5>(gy (36) 

i=0 

"2 

d(a,g) = i;6 i (gy (37) 

where the complex coefficients c«(g) , are affine functions of the uncertain parameters 

q = (qx , q 2 , , qif , respectively. 

In the sequel, we will suppose that D c is a connected set. Note that Hurwitz stability 
and Schur stability are special cases of D-stability. 

Lemma 5 

For any fixed complex numbers Zqi , Zq 2 & C , suppose the polynomial family {zqi7i(s, q) + 
z 02 d(s, q) | q G Q} has a fixed order. Then 

{z 01 n(s, q) + z 02 d(s, q) \ q E Q} C D 

(38) 

{z Q1 n(s, q) + z 02 d(s, q) \ q e Q E } C D 

Proof: Since the coefficients of z in(s, q) + z 02 d(s, q) are also affine functions of q = 
(qx , q 2 , , qi) T , the result follows directly from the Edge Theorem[|, 

For notational simplicity, define 
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g{s,q) = a m [n(s,q)] m + a m -i[n(s, q)] m 1 d(s,q) 

+a m ^[n(s,q)] m ~ 2 [d(s,q)] 2 + 

+a 2 [n(s, q)] 2 [d(s, q)] m ~ 2 + ai n(s, q)[d(s, q)}^ 1 

+a [d(s,q)] m 

where a& e C , k — 0, 1, , m . 

Theorem 3 

Consider the polynomial family 

{g(s,q) \qeQ} 

Suppose it has a fixed order. Then 

{g(s,q) | q e Q} C D 

{g(s,q) \ q eQ E }cD 

Proof: Consider the polynomial 

q(z) = a m z m + a m _iz m_1 + a m _ 2 ^ m ~ 2 + + a-2Z 2 + a\z + a 

Let r = max{/c | afc ^ 0}. Then can be expressed as 

q(z) = a r (z - Zi)(z — z 2 ) (z — z r -x)(z - z r ) 

where z± , z 2 , , 2 r - 1 , z r G C. Hence, we have 
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{fl(s> <?) I 9 6 Q} C D g(s,q) £ -D, Vq S Q 

<)P {-m " + m " + • • • 

a r [d(s, q)]™~ r [n(s, q) — z\d(s, g)][n(s, q) — z 2 d(s, q)] ■ ■ ■ 

■ ■ • [n(s, q) — z r _ 1 d(s, q)][n(s, q) ~ z r d(s, q)] £ D . Vq G Q 



n(s,q) - z k d(s,q) G D , 

k = 1,2, ,r - l,r, Vg G i 

- z k d(s,q) G D , 

k = 1,2, , r - l,r , Vg G i 

and <Z(s, g) G D , Vg G Q 



n(s,g) - z k d(s,q) G D , 

fe = 1, 2, ,r - l,r, Vg G Q E 

n(s, g) — z k d(s, g) G £> , 

fc=l,2, ,r-l,r, VgGQ B 

and d(s,q) G D , Vg G Q B 



r < m 



r < m 



(44) 



a r [d(s, q)] m r [n(s, q) — 2id(s, g)][n(s, q) — s 2 (i(s, q)] • ■ ■ 

■ ■ ■ [n(s, q) — s r _id(s, q)][n(s, q) — 2 r d(s, q)] £ D , \/q £ Q E 

W'.«)l m {-r[«-*l] 

»)] m {-«. ™ + -m-i [^fffy] m " + ' • • 



g(s, q) G D , Vq S Q E 



Vq 6 Q E 
{S(s, q) | q € Q E } C D 



This completes the proof. 

Remark. Theorem 3 reveals that, for a class of polynomial family with nonlinearly corre- 
lated perturbations, D-stability of the entire family can be ascertained by only checking 
one-dimensional edge polynomials in this family. 
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(45) 



(46) 



3.2 Extension to Polynomial Matrix Families 

Consider the uncertain polynomial matrix 

" 2n(s, q) + 3d(s, 9) 3n(s, q) + 4d(s, q) 

M(s,q)= 4n(s , q) + 5d(s , q) 2n(s, q) 

9d(s, q) 6n(s, q) 5n(s, q) + Gd(s, q) 

it is easy to see that 

det[M(s, q)} = 16[n(s, q)] 3 + 176[n(s, g)] 2 d(s, g) 
+279n(s, g)[d(s, g)] 2 + 90[d(s, g)] 3 

By Theorem 3, we have 

{det[M(s, q)} \ q G Q} C D 

(47) 

{det[M(s,g)] |gGQ £ }cZ) 

Namely, robust D-stability of the entire polynomial matrix family can be ascertained by 
only checking one-dimensional edges. More generally, for any uncertain polynomial matrix 
of the form 

M(s, q) = [a i:j n(s, q) + /%d(s, q)] nxn (48) 

it is easy to see that the above edge result also holds. Moreover, if n(s, q) , d(s, q) are 
replaced by interval polynomial families Y u , r„ or A u , A„ as defined in the last section, 
then Kharitonov-like results can be established for robust Hurwitz stability of the corre- 
sponding polynomial matrix families. 

Theorem 4 

Consider the polynomial matrix family 

M(5 u (s), 8 v (s)) = hj$u(s) + r]iA(s)}nxn (49) 
where S u (s) G A u , S v (s) G A„, and jij,T]ij,i,j = 1,2, , n are complex numbers. Then 

{det[M(5 u (s), 5 v (s))\ \ 5 u (s) G A M , 5 v (s) G A„} c H 

{det[M(Kt u (s) , K^(s))] | i,j = 1,2,3, 4}U (50) 
{det[M(K- u (s) , Kr^s))] \ i,j = 1,2,3,4} C H 
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4 Some Applications 



A proper transfer function 44 is said to be strictly positive real, denoted by 44 G SPR, 
if 



1) q(s) G H 

2) 3?j|g > , Vcg> g i? 

Suppose have positive leading coefficients. Then, it is easy to see that 

p(s) 



G 5Pi? Ap 2 (s) + (1 - A)<f (s) G H, X G [0, 1] 



Now consider the proper interval transfer function family 



In order to have 



T 



Pu{S) 

Pv(s) 



Pu{S) 
Pv{s) 



Pu{s) e r u , p v (s) g r\ 



g spr , p u (s) g r u , p v (s) g r„ 



we must have 



(51) 



(52) 



(53) 



(54) 



\TI + (1-\)T 2 V CH,\E[0,1} (55) 

Since Xz 2 + (1 — A) has purely imaginary roots. By Theorem 1, we only need to have||, 
|J,1 







- A) 


[K v 4 (s)} 2 EH,Xe 


[0,1] 


(56) 


X[K%[ 




-A) 


[K v 3 (s)] 2 eH,Xe 


[0,1] 


(57) 


X[K-( 




-A) 


[K^s)} 2 eH,Xe 


[0,1] 


(58) 


X[K%{ 


*)] 2 + (l 


-A) 


[K v 2 (s)} 2 eH,Xe 


[0,1] 


(59) 


\[K?{ 


*)] 2 + (l 


-A) 


[K v 3 (s)] 2 eH,Xe 


[0,1] 


(60) 


X[K^[ 


*)] 2 + (l 


-A) 


[Kl{s)} 2 eH,Xe 


[0,1] 


(61) 


X[K-( 


*)] 2 + (l 


-A) 


[KV(s)} 2 eH,Xe 


[0,1] 


(62) 


X[K U 4 ( 


*)] 2 + (l 


-A) 


[K^s)} 2 eH,Xe 


[0,1] 


(63) 



11 



Equivalently 

*?(«) jgto ^4(«) 

g£W g£W *W gj&) c Cpn 

Namely, in order to guarantee that every member of the interval transfer function family 
T is strictly positive real, we only need to check eight specially selected vertex transfer 
functions. That is 



Puis! 

Pv(s) 



e SPR , \J Pu (s) e T u , V Pv {s) e T v 

(65) 



K?(s) K%{s) KX(b) K%(s) 
Kl(s) ' iq(e) ' K-(s) ' K$(s) 



g£(f) K 2(s) Kg (s) K^s) qpR 
K$(s) ' Kl(s) ' K%(s) ' KZ(s) C orlx 



which is consistent with the result of Chapellat et al|J, and WangJTU 
Moreover, for any 7 G R , in order to have 

7 + ^4 e spr > pm e r - > Pv( s ) e r - ( 6e ) 

we must have 



A[ 7 r, + r u ] 2 + (1 - X)T 2 V cH,\e [0, 1] (67) 

Since A (7 + z) 2 + (1 — A) has roots either at first and fourth quadrants (when 7 < 0) or at 
second and third quadrants (when 7 > 0). By Theorem 1, we only need to check twelve 
vertices to guarantee robust stability 0, ||, [7], ||. Namely, in order to guarantee that 

7 + ^44 e spr > pm g r « > g r - ( 6g ) 

Pv[S) 

we only need to check the same property for twelve specially selected vertex transfer 
functions. 

Remark. The above result can be easily extended to the case of complex interval transfer 
function family. Namely, every member in the complex interval transfer function family 
is strictly positive real, if and only if, sixteen specially selected vertex transfer functions 
in this family are strictly positive real [HI . 



5 Robust Sensitivity Functions 

Denote the m-th, n-th (m < n) order real interval polynomial families K g (s), Kf(s) as 

m 

K g {s) = {g{s)\g{s) = e ikM* = 0,1, M, (69) 

i=0 
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K f( s ) = {f( s )\f( s ) = ^2aiS l ,ai G [ai,ai\,i = 0, 1, ,n}. (70) 

i=0 

For any f(s) G Kf(s) , it can be expressed as 

f(s) = a f (s 2 ) + s(3 f (s 2 ), (71) 

where 

«/(s 2 ) = a + a 2 s 2 + a 4 s 4 + a 6 s 6 + , (72) 

(3 f {s 2 ) = ai + a 3 s 2 + a 5 s 4 + a 7 s 6 + (73) 

Obviously, for any fixed to G R, a.f(—L) 2 ) and ujf3f(—u} 2 ) are the real and imaginary parts 
of f{jw) G C respectively. 

For the interval polynomial family Kf(s), define 

ci!^(s 2 ) = ap_ + a^s 2 + CI4S 4 + a^s 6 + , (74) 

af\s 2 ) = oq + a2_s 2 + ajs 4 + OqS 6 + , (75) 

pf\s 2 ) = ai + a^s 2 + a 1 s 4 + a^s 6 + , (76) 

pf\s 2 ) = (H + a A s 2 + aEs 4 + a 1 s 6 + , (77) 

and denote the four Kharitonov vertex polynomials of Kf(s) as 

f lJ (s) = af(s 2 ) + s(3f(s 2 ), i,j = 1,2 (78) 

For the interval polynomial family K g (s), the corresponding a^'(s), /3jf>(s) and gij(s) G 
K g (s) can be defined analogously. 



Lemma 6[13| 



For any fixed u G R, f(s) G Kf(s), we have 



a 



{1 \-co 2 ) < a f {-u 2 ) < afi-u 2 ), (79) 



Pf\-J)<Pf{-J) (80) 

Lemma 7||TT| (Zero Exclusion Principle) 
For the n-th order polynomial family 

f(s,T)=:{f(s,t)\tET}, (81) 

where T is a bounded connected closed set, and the coefficients of f(s,t) are continuous 
functions of t, then f(s,T) G H if and only if 
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1) there exists t* G T, such that f(s,t*) G H; 

2) 0$f(ju,T),VueR. 

Consider the strictly proper open-loop transfer function 

and suppose the closed-loop system is stable under negative unity feedback. Denote its 
sensitivity function as 

„ i /w (83) 



i+p /(»)+»(») 

Apparently, we have 



|S1|oo>1 (84) 



For notational simplicity, define 



Jiijii2j 2 (s) = 9i 1 ji( 8 ) + i 1 + Se3<> )fi2j2(s), S £ (0,1), ii,ji,i 2 ,j2 = 1,2, 9 G [— 7T,7T]. 

(85) 

Lemma 8 

Suppose <7(s) + /(s) G if. Then, for any 7 > 1, we have 

Halloo < 7^ g(s) + (l + -e jd )f(s)EH, We[-n,n]. (86) 

7 

Proof: Necessity: Since + /(s) G H and II /(J) +g ( s ) lloo < 1, by Rouche's Theorem, we 



know that 



W) + f(s)} + -e' e f(s) G H, W9e [-tt, tt] (87) 

7 

10 r . 



Sufficiency: Now suppose on the contrary that H^Hoo > 7, namely, 1 1 j^y^ 
Since | f^)+g( s ) \ s=ju\ * s a contiunous function of to, and since 



, , 1 



lim I 7 | -J = - < 1 (88) 

w - >0 ° f(s)+g(s) 7 



there must exist u;q such that 



7 /(*) 



7W + ffW 

Therefore, there exists 9 G [— 7r, tt] such that 



l*=J^0 I 



(89) 
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(90) 



Lemma 9 

For any 5 G (0, 1), 9 G [—it, 7r], we have 



W(s) =: {(/(a) + (1 + G /(a) G C H 



(91) 



•fllll) ^1212, ^2222, ^2121) ^1112) <fl222) ^2221) ^2 111) ^1211) ^2212) ^2122) ^1121 G if (92) 



Proof: Necessity is obvious. To prove sufficiency, note that W(s) is a set of polynomials 
with complex coefficients, and with constant order n. By Lemma 7, it suffices to show 
that 



where dW(ju) stands for the boundary of W(ju) in the complex plane. 

To construct dW(ju>), note that arg(l + 5e J ) G (— | , |). Suppose now uj > and 
arg(l + <5e je ) G [0, |). Then by Lemma 6, we know that K g (ju), Kf(ju>) are rect- 
angles with edges parallel to the coordinate axes. The four vertices of K g (juj) are 
gii(juj),g 1 2(juj),g2i(juj),g22(ju), respectively; and the four vertices of K f (ju) are fu(juj), 
fuiju), f2i(jv), /22CM, respectively. (1 + 8e> e )Kf(ju) is generated by rotating K f (Ju) 
by arg(l + 5e^ e ) counterclockwisely, and then scaling by |1 + 8e J |. Thus, W(juj) = 
Kg(ju) + (1 + 5e^ d )K f{ju) is a convex polygon with eight edges. These edges are parallel 
to either the edges of K g {juj) or the edges of (1 + SeP )Kf(Ju). Therefore, their orienta- 
tions are fixed (independent of a;). The eight vertices of W(ju) are (clockwisely) Jun(ju), 
JiiuUuj), Ji2i2(j^), J1222UU), J 2 222(j^), ^222i (juj), J 2 i2i(j^), ^2in (ju), respectively. 

Now suppose on the contrary that there exists u > such that 



g" W(juj), MujeR 
Since g" W(j Uoo) for sufficiently large u^, we only need to show that 



(93) 



OgdW(ju), WueR 



(94) 



G dW(ju; ) 



(95) 



Without loss of generality, suppose 



G {AJmiCM,) + (1 - A) J m2 0'wo)|A G [0, 1]} 



(96) 



Namely, there exists A G (0, 1) such that 



A Jiiii(j^o) + (1 - A ) Jm 2 (M)) = 



(97) 



Since Jim(s), Juuis) G H, we have 
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Thus [HI 



d d 
— arg Jim (jo;) > 0, — axg J m2 (ja;) > 
duo duo 



(9* 



arg[Jiii 2 (ju;) - Jmi(j^) 



Ao)^j arg Jim (jw) | 



d 



u=u + Ao^arg J 1 ii 2 (^)| w =w > 



(99) 
(100) 



This contradicts the fact that the edges of W(juo) have fixed orientations. Thus 



£ dW(juo) (101) 

Suppose now uo < and arg(l + ^e- 761 ) G (— f , 0]. Then K g (juo), (1 + 5e^ e )Kf{juo) are the 
mirror images (with respect to the real axis) of the corresponding sets in the case of uo > 
and arg(l + Se^ 9 ) G [0, |p. Therefore, following an identical line of arguments, we have 



£ dW(juo) (102) 

The cases when uo > and arg(l+5e je ) G (— |, 0] and when uo < and arg(H-5e- ?e ') G [0, |) 
are also symmetric with respect to the real axis. Hence, we only need to consider the for- 
mer case. In this case, K g (juo), Kf(juo) are rectangles with edges parallel to the coordinate 
axes. (1 + 8e? )Kf{juS) is generated by rotating Kf(juo) by | arg(l + <5e j6, )| clockwisely, and 
then scaling by |1 + 5e j6 '|. Thus, W(juo) = K g (ju>) + (1 + 5e^ d )Kf(juo) is a convex polygon 
with eight edges. These edges are parallel to either the edges of K g (ju) or the edges of 
(1 + 5ei e )Kf(juo). Therefore, their orientations are fixed (independent of uo). The eight 
vertices of W(juo) are (clockwisely) J lul (joo), JmiO'w), J^uU^), J2212UU), J 2 222(j^), 
J2i22(j^), J2i2i(j^)) Jii2i(j^), respectively. Thus, following a similar argument, we have 



£ dW{juo) (103) 

This completes the proof. 

The following theorem shows that, for an interval system, the maximal norm of 
its sensitivity function is achieved at twelve (out of sixteen) Kharitonov vertices. 

Theorem 5 

Suppose gij(s) + fij(s) G H, i,j = 1,2. Then 

max{ll 7(^T^) l|oo|g(s) G K ^'f^ e K ^ = ( 104 ) 

max{|| J UKW2J2) = (1111), (1212), (105) 



(2222), (2121), (1112), (1222), (2221), (2111), (1211), (2212), (2122), (1121)} (106) 
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Proof: Since g%j{s) + fij(s) G H, i,j = 1,2, by Kharitonov's Theorem[|IJ, we know that 
K g (s) + K f (s) C H. Let 

7l = max{|l /(a) J |ocb(a) e K g (s)J(s) G K f (s)} (107) 



700 + sOO 1 



72 = max{H - . J |oo|(W2j 2 ) = (1111), (1212), (108) 



(2222), (2121), (1112), (1222), (2221), (2111), (1211), (2212), (2122), (1121)} (109) 
Then apparently 

7i>72>l (HO) 

Now suppose 71 ^ 72, namely, 71 > 72 . Then there exists 70 such that 71 > 
70 > 72- Thus, for any (w 2 J 2 ) G {(1111), (1212), (2222), (2121), (1112), (1222), (2221), 
(2111), (1211), (2212), (2122), (1121)}, we have 



<7o (HI) 



Hence, by Lemma 8, we have 



(^(s) + (1 + -e^)4 J2 ( S ) € H, \/6 G [-7r,7r] (112) 
7o 

By Lemma 9, we know that 

{g(s) + (1 + le3«)/( a )|^( s ) e K g (s)J(s) G K f (s)} C H, W G [-vr, vr] (113) 
7o 

Therefore, by Lemma 8, for any g(s) G K g (s), f(s) G Kf(s), we have 

M /00 



< 7o (H4) 



Namely 



max{|| - /(s) r |[ocb(g) G K g (s),f(s) G < 70 (115) 



That is, 71 < 70, which contradicts 71 > 70 > 72- This completes the proof. 
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6 Conclusions 



Some Kharitonov-like robust Hurwitz stability criteria have been established for a class 
of complex polynomial families with nonlinearly correlated perturbations. These results 
have been extended to the polynomial matrix case and non-interval D-stability case. 
Applications of these results in testing of robust strict positive realness of real and complex 
interval transfer function families have also been presented. 
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